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Abstract 



Gravitational wave detectors in the mHz band (such as LISA) will observe thousands compact binaries in the galaxy which can be used to better understand the structure of the Milky Way. To test the effectiveness of 
LISA to measure the distribution of the Galaxy, we simulate the Close White Dwarf Binary (CWDB) gravitational wave sky using different models for the Milky Way. To do so, we have developed a Galaxy density 
distribution modeling code based on the Markov Chain Monte Carlo method. The code uses different distributions to construct realizations of the Galaxy. We then use the Fisher Information Matrix to estimate the 
variance and covariance of the recovered parameters for each detected CWDB. This is the first step towards characterizing the capabilities of space-based gravitational wave detectors to constrain models for Galactic 
structure, such as the size and orientation of the bar in the center of the Milky Way. 


Introduction 


The Disk’s Galactic Distribution 


One of the science goals of LISA is to understand and map the structure of the Galaxy and the distribu- 
tion of close white dwarf binaries (CWDBs). One way of achieving this goal is to obtain sky locations 
for individually resolvable CWDBs along with distances for chirping CWDBs. The Mock LISA Data 
Challenges (MLDCs) have shown that we can expect to resolve ~ 20, 000 binaries. A large number of 
these systems (those at frequencies > 4 mHz) will be distributed homogeneously throughout the galaxy. 
In order to estimate the degree to which LISA observations can constrain Galaxy structure models, we 
anticipate modeling the errors expected from typical data analyses and applying the results to a variety 
of Galaxy models spanning astrophysically reasonable structure parameters. We model the CWDB pa- 
rameter estimation errors from different Galaxy models containing a bar and from this recovered catalog 
of CWDBs determine if LISA is able to detect this Galactic structure. 



Figure 1: LISA detector (Right) and a White Dwarf binary star propagating gravitational waves(Left) 

Process 


We begin with a catalog of approximately 27 million CWDBs containing derived from population syn- 
thesis simulations of the Galaxy. Using this distribution we can Monte Carlo over each binaries position 
in the galaxy, and orientation with respect to the LISA detector, to produce different Galaxy models ac- 
cording to the different density profiles shown in Eq 1 - 8. A Markov Chain Monte Carlo algorithm was 
used to generate new instances of each Galaxy model - in other words, with every run of the code we 
generate a new Galaxy. 

After obtaining different Galaxy models we determine which binaries are detectable using a signal-to- 
noise (SNR) threshold of 7. From this subset of detectable binaries, the Fisher Information Matrix is 
used to estimate each binary’s parameter uncertainties. A further selection is made on those which are 
resolved by LISA to within 1 deg 2 on the sky, and with there luminosity distance d i measured to within 
10%. This subset of spatially-resolved binaries will be used in future work to constrain the distribution 
of the Galaxy. 


Galaxy using bar nodel El uith 1 nillion sanples Galaxy using bar nodel El uith 1 nillion sanples 
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We start to populate the disk of the Galaxy by selecting R and z from the probability density function 
(PDF) 

p(R, z ) = pod e~^^ d e~ z ! Zd (1) 

with scaleheight zy, scalelength R ( j and normalization constant po jj . 

The galactic disk lies in the x-y plane. The distribution of stars is assumed to be cylindrically symmetric, 
so the azimuthal angle 9 is uniformly drawn from 0 < 9 < 2tt. We then compute the cartesian galactic 
coordinates of each binary via x — R cos(9), y — Rsm(9), and R — + p 2 . 

The Milky Ways disc is usually considered to have two major components: a thin disc and a thick disc 
(e.g. Gilmore & Reid 1983). As described in Juric et al.(2008), the scalelengths R ( j for the thin and thick 
discs are 2.6 kpc and 3.6 kpc, respectively, with a quoted uncertainty of 20 per cent in each case, and the 
best-fitting values for the scalengths are zy+j = 300 pc and zy thick = 900 pc. In this treatment, we only 
consider a single disk distribution with scale height zy = 300 pc. 

The Bulge and Bar’s Galactic Distribution 

We populate the Bar of the Galaxy by using the probability density functions used in the analytic models 
described by Dwek et al. (1995). The density functions are the following: 
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and the normalization for the bulge’s distribution function is 
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The coordinate system has the origin at the GC, with the xy plane defining the mid-plane of the Galaxy 
and the z-direction parallel to the direction of the Galactic poles. 


Model 
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g 2 

-7.68 

24.49 
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473.28 

359.91 

Ei 

1.95 

23.82 

1901.43 

626.89 

324.93 

E 2 

0.97 

26.43 

986.60 

356.65 

260.88 

Pi 

2.19 

25.30 

1810.98 

609.40 

478.93 

P2 

3.67 

25.16 

3513.34 

1160.72 

927.54 


Table 1 : Best-fitting parameter values for all density models (equations 2 tolO) with Bar orientations (/? and a ) in degrees 
and scale lengths (:ro, yo and zq) in parsecs. 
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Figure 2: Front view of a Galaxy model using best bar model El (top left) and the bulge (bottom left) with their respective top 

views (top right and bottom right respectively). 


References 


5. Conclusion 

From the different Galaxy models used, it was observed that many of the LISA resolvable CWDBs will be 
coming from the bar and bulge as seen in figure 3. This due to the high concentration of sources located 
in these Galactic structures. From this results it is intended in future work to constrain the distribution 
of the Galaxy. Other Galactic structures can be implemented in the models such as the addition of spiral 
arms and globular cluster for more realistic models (as many papers suggest). 

LISA catalog of CHDBs using nodel El LISA catalog of CHDBs using the Bulge nodel 
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